The Collinear Magnetic Phases of the Geometrically-Frustrated Antiferromagnet 

CuFe02: The Importance of Stacking 
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The correct stacking of hexagonal layers is used to obtain accurate estimates for the exchange and 
anisotropy parameters of the geometrically- frustrated antiferromagnet CuFe02. Those parameters 
are highly constrained by the stability of a collinear metamagnetic phase between fields of 13.5 and 
20 T. Constrained fits of the spin-wave frequencies of the collinear ttii phase below 7 T are used 
to identify the magnetic unit cell of the metamagnetic tttii phase, which contains two hexagonal 
layers and 10 Fe^^ spins. 
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PACS numbers: 75.30.Ds, 75.50.Ee, 61.05.fg 



Because of their rich magnetic phase diagrams, 
geometrically-frustrated antiferromagnets have long oc- 
cupied an important place in condensed-matter physics 
[l| . The antiferromagnetic interactions between the Fe'^+ 
spins of CuFc02 are geometrically frustrated within each 
hexagonal plane since no spin configuration can simul- 
taneously minimize the coupling energies of all three 
neighbors around an equilateral triangle. Unlike for 
many geometrically-frustrated antiferromagnets, quan- 
tum fluctuations about the magnetic ground states of 
CuFe02 can be safely neglected due to the large S = 5/2 
spins. Whereas geometric frustration often leads to mag- 
netic phases with non-coUinear spins and complex unit 
cells, magnetic anisotropy perpendicular to the hexag- 
onal planes in CuFe02 produces two different collinear 
magnetic phases. The ttii phase [1, 0] sketched in 
Fig. 1(a) is stable up to the field Bd ~ 7 T . Between 
Bc2 ~ 13.5 T and Scs ~ 20 T, another collinear phase 
with a net moment of 1 /iB per Fe'^"'" ion 0, Q has been 
assumed to resemble the tttii phase shown in Fig. 2 for 
type B stacking, with 5 Fe^+ spins per unit cell. In- 
commensurate and non-collinear phases were identified 
between Bd and Bc2 and above -Bc3 0, la] ■ 

Previous efforts to understand the collinear magnetic 
phases [1, [j, la] and to estimate the exchange and 
anisotropy parameters I7l| of CuFe02 made the simpli- 
fying assumption that the hexagonal layers were stacked 
sequentially on top of each other. We now demonstrate 
that an accurate determination of the Heisenberg param- 
eters must employ the correct stacking of the hexagonal 
layers. We also show that the stability of a metamagnetic 
phase between Bc2 and Bci, [1, M, strongly constrains 
those parameters. Whereas earlier work [3] assuming a 
sequential stacking was unable to explain the observed 
spin-wave (SW) frequencies of the zero-field twins, re- 
alistic magnetic stackings are now used to explain all 
features of the low-field collinear phase and to identify 
the magnetic unit cell of the high-field collinear phase in 
CuFe02. 

The observation of collinear magnetic phases that are 



fully polarized along the ±z directions at low tempera- 
tures led to the assumption [3,[g] that the Fe'^^ spins were 
"Ising-like." However, measurements of the zero-field SW 
frequencies y, ^ plotted in Fig. 1(b) reveal SW gaps of 
only about 0.9 meV at wavevectors [H, H,L — 3/2) with 
H = 0.21 and 0.29, on either side of the ordering wavevec- 
tor Q — (1/4,1/4,3/2). If the spins were truly "Ising- 
like," then the SW frequencies would be much higher and 
they would not exhibit a significant dispersion along the 
(0, 0, L) direction [2, |9l] perpendicular to the hexagonal 
planes. With little change in wavevectors, the SW gaps 
are reduced either by an applied field along the z axis or 
by the substitution of nonmagnetic AP+ ions for Fe^+. 
Above the field Bd \% l5| or an Al concentration of about 
1.6% Q, the SW gaps vanish, the magnetic ground state 
becomes non-collinear, and the crystals display multifer- 
roic behavior 0, [HI, [ll ■ 

Assuming that the hexagonal planes stack sequentially, 
we recently fit Q the SW frequencies of pure CuFe02 to 
the predictions of the Heisenberg model 
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which includes single-ion anisotropy D and a magnetic 
field B. For "Ising-like" spins, D would be much greater 
than the exchange parameters Jij. In a further simpli- 
fication, we ignored the very small (< 0.4%) distortion 
of the hexagonal plane p, [lj| below the Neel temper- 
ature. While this distortion breaks the symmetry be- 
tween the {H,H,0), (iJ, 0,0) and (0, i/,0) directions, 
thereby favoring the ttii phase with wavevector Q over 
its twins, it can produce only a very small change in the 
exchange parameters and hence in the SW frequencies. 
Despite these simplifications, the SW dispersions evalu- 
ated along the (H, H, 3/2) and (0, 0, L) axis agree quite 
well with inelastic neutron-scattering measurements [t!]. 
However, we were unable to fit the frequency of the two 
twins with wavevectors rotated ±7r/3 away from Q in 
the {H, K, 3/2) plane. Without attempting to fit the 
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FIG. 1: (a) The low-field spin configuration (up spins are 
empty and down spins are filled circles) in each hexagonal 
plane, with the four inequivalent spins a, b, c, and d. Both the 
n (sofid) and n+1 (dashed) layer are shown with the exchange 
parameters indicated, (b) The fit of the S W frequencies along 
the {H, H, 3/2) axis using the exchange and anisotropy pa- 
rameters given in line in of Table I. Open squares give the 
frequencies of the twins and solid circles the frequencies of the 
main SW branch with ordering wavevector at H = 1/4. 



twins, we obtained the exchange and anisotropy param- 
eters given in line i of Table I, where Jpm or J^m are the 
mth nearest-neighbor exchange parameters within each 
hexagonal plane or between adjacent planes. 

To better understand the metamagnetic phase between 
Bc2 and Bcs, we have recalculated the SW frequencies of 
the ttii phase below Bd using the realistic magnetic 
stacking of the hexagonal layers shown in Fig. 1(a). All 
other stackings of the ttii layers have higher coupling 
energies. Because spins a, 5, c, or d experience the same 
local environment on every layer, the magnetic unit cell 
still contains only 4 sublattices (SLs). The first few ex- 
change pathways Jpm and Jzm are indicated in Fig. 1(a). 

The SW frequencies are evaluated using a Holstein- 
Primakoff (HP) l/S expansion about the classical limit. 



On the spin-up a and b sites, we replace Siz 



alcti, 



TABLE I: Heisenberg parameters of CuFe02 obtained from 
fits of the zero-field SW frequencies. Line i assumes sequential 
stacking of the hexagonal layers [31 , ii and iii use the realistic 
stacking in Fig. 1(a) while iii also constrains the parameters to 
stabilize the coUinear phase between Bc2 and Bcs- Exchange 
and anisotropy parameters are in meV; T^^^ is in K. 
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FIG. 2: Three types of magnetic stacking that satisfy the 
conditions for local stability of the metamagnetic phase. In 
type A stacking, the stacking patterns on the left and right 
alternate. 



(cki — ai or hi). On the the spin-down c and d sites, we 
replace S^z = -S + -fjJ^, S+ = V2Sj}, and 5" = V2Sj^ 
{li = Ci or di). The SW frequencies Wk at wavevector k 
are then obtained by solving the equations-of-motion for 
the vectors Vk = (ok, &k,cj^,rfk) and vj^. The equation- 
of-motion for Vk niay be written in terms of the 4x4 
matrix M(k) as idv^/dt = ~[H, Vk] = M(k)vk with SW 
frequencies given by the condition Det( M (k) — muD = 0. 
Only positive frequencies Wk > are retained. 

As expected for a collinear antiferromagnet and shown 
schematically for any wavevector in Fig. 3(a), each of the 
SW branches is linearly split by a magnetic field. The 
lowest SW frequency with wavevector (0.21,0.21, 1.5) or 
(0.29,0.29,1.5) wiU vanish at the field 0.9 meV/2fiB ~ 
7.7 T, which is slightly larger than Bd- With the correct 
stacking of the hexagonal layers, the parameters in line 
ii of Table I are obtained by fitting the SW frequencies 
of the main branches along the {H, H, 3/2) and (0, 0, L) 
axis as well as the SW frequencies of the twins evaluated 
along (iJ,iJ, 3/2). 

There are several possible collinear metamagnetic 
phases with a net moment of 1 ^b per Fe^^ ion and with 



elastic peaks at wavevectors {m/5,m/5,0) in the L = 
basal plane [j|. Two configurations are possible in each 
hexagonal plane: the TTTii pattern sketched in the lower 
left of Fig. 2 and the ttiti pattern sketched in the lower 
right. Depending on the stacking, the magnetic unit cell 
of the metamagnetic phase may contain either 5 or 10 
magnetic ions. For example, type A stacking of tttii 
layers in Fig. 2 contains 10 SLs, while type B stacking of 
TTTii layers and type C stacking of TTiTi layers contain 
5 SLs. In type A stacking, the local environments of spin 
a on layer n and spin a' on layer n + 1 are different: a 
is coupled by J^i to three up spins on layer n + I while 
a' is coupled by Jzi to two up spins and one down spin 
on layer n + 2. In types B and C stacking, the spins 
in layer n ~\- 1 are obtained from those in layer n by the 
displacement — V3y/3. For a 5 or 10 SL stacking, the 
matrix M(k) that enters the equations-of-motion for the 
SW frequencies is 5 or 10 dimensional and the 5 or 10 
SW branches s must be solved numerically for every k. 

Two conditions must be satisfied for the local stabil- 
ity of a metamagnetic phase. First, the SW frequencies 
Wj^ must all be real. This condition is independent of 
the magnetic field, which only shifts the frequencies by 
zt2fj,BB, and is not always satisfied because M(k) is not 

Hermitian. Second, the SW weights W^^ that appear as 
coefficients of the delta functions in the spin-spin corre- 
lation function 

Sik,u;) = 1 /die— '5]e^''-(^^-^-){(5+5-(i)) 

+ {srs+it))} = J2 <'^^i^ - 4"^) (2) 

s 

must all be positive. Those weights are most easily evalu- 
ated by expanding 5(k, w) within the HP formalism and 
then solving the equations-of-motion for the spin Green's 
functions. An equivalent but much easier way to guaran- 

(s) 

tee that the weights W^ are positive is to examine the 
field-dependence of the SW frequencies. For a stable 5 
SL coUinear phase, 3 of the 5 SW modes must linearly in- 
crease with field while 2 must linearly decrease, as shown 
in Fig. 3(a). For a stable 10 SL coUinear phase, 6 of the 
10 SW modes must linearly increase and 4 must linearly 
decrease with field. If this condition is violated for any 

(s) 

k, then some of the weights W^^ will be negative and 
the phase will be unstable. 

Unfortunately, the exchange and anisotropy parame- 
ters given by lines i and ii of Table I do not satisfy both 
conditions for the local stability of any possible stacking 
of TTTii or TTiTi layers between the fields Bc2 and Bcs- 
In other words, fits to the SW frequencies of the zero- 
field TTii phase are inconsistent with the existence of a 
coUinear metamagnetic phase. 

This inconsistency may be eliminated by fitting the 
zero-field SW frequencies of the TTii phase while simul- 
taneously constraining the exchange and anisotropy pa- 
rameters to stabilize a metamagnetic phase between Bc2 
and Bcs- We emphasize that this constraint utilizes only 
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FIG. 3: (a) Schematic field dependence of the SW frequencies 
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j^ in 4 or 5 SL phases, (b) The predicted elastic intensities 
\F{m)\'^ normalized by the Fe'^"'" form factor f{m)^ versus m 
(H = m/5 along the {H,H,0) axis) for stackings A, B, and 
C of the metamagnetic phases. Inset are the experimental, 
normalized intensities versus m. 



the observed stability of the metamagnetic phase over a 
range of magnetic fields and not the measured SW fre- 
quencies of that phase. The three phases shown in Fig. 2 
are the only ones that satisfy both conditions for local 
stability when the exchange and anisotropy parameters 
are obtained from constrained zero-field fits of the SW 
frequencies. 

To determine which of these three phases is observed, 
we evaluate the magnetic structure factor F{m) for the 
elastic peaks in the L = basal plane at wavevectors 
{H,H,0) with iJ = m/5: 



F^{m) = f{m)e 



2TriH 



(3) 



F^{m) = /(m)e2'^*^{ 



2i sin 27rH + 2i sin inH + 1 



2 cos 2t:H - 2i sin 4t:H - 1 



}, (4) 
}, (5) 



where /(to) is the magnetic form factor of each Fe'^^ ion. 
Notice that F{Q) = /(O) for all three possible phases. 
The normalized intensities \F{m)\'^ / /{m)'^ are plotted 
versus to, in Fig. 3(b). When the magnetic moments of 
6 adjacent layers are summed, stacking A produces the 
pattern OTOOO along x so that \F'^{m)\'^/f{mf = 1 is 
constant. If the stacking pattern on the left or right top 
panel of Fig. 2 were continued indefinitely rather than al- 
ternating, then the resulting phase would have no elastic 
peaks in the L = basal plane. The layer sums of stack- 
ings B or C produce TTTii or TTiTi patterns along x, 
causing \F{m)\'^ / /{m)'^ to change by a factor of 10.5 as 
m increases from 1 to 2. 

For comparison, the experimental results 10] for the 
elastic intensities are plotted in the inset to Fig. 3(b). The 
normalized intensity \F(m)\'^ / f{m)'^ for to = 1 through 
4 is constant to within about 1%. Therefore, only type A 
stacking of TTTii layers with a 10 SL unit cell is possible. 

The exchange and anisotropy parameters associated 
with stacking A are given on line Hi of Table I. As in our 



original fits [?!], \Jp3\ > \Jp2\ but Jp4 is negligible. Since 
Jz3 is comparable to J^i, even longer-ranged interactions 
between neighboring planes might exist. All of the inter- 
actions Jzm between adjacent planes are much smaller in 
magnitude than the interactions Jpm [m < 4) within a 
plane. Using these parameters, the fits of the main and 
twin SW branches are plotted along the {H,H,3/2) axis 
in Fig. 1(b). 

Constraining the fits of the zero-field SW frequencies 
to produce a stable metamagnetic phase has a substan- 
tial effect on the exchange and anisotropy parameters. 
For example, Jpi is reduced by about 70% from line ii 
to line in of Table I. While a wide range of parameters 
can provide reasonable fits to the zero-field SW data, de- 
manding that a metamagnetic phase is stabilized between 
Bc2 and Bcs considerably narrows the possible range of 
those parameters. Also notice that the mean-field tran- 
sition temperature Tj!f^ listed in line in of Table I is 
much closer to the measured transition temperature of 
14 K [J| between partially-disordered and paramagnetic 
phases than the transition temperatures of the uncon- 
strained fits in lines i and ii. Of all three fits, line Hi 
produces a crystal-field environment that is most "Ising- 
like," with the anisotropy D about the same size as the 
nearest- neighbor exchange Jpi. The difference between 
the parameters in lines i, ii, and Hi of Table I underscores 
the danger of using even an extensive set of SW measure- 
ments for a single magnetic phase to fix the parameters 
of a Heisenberg model. 

Surprisingly, the high-field coUinear phase is the 10 SL 



phase sketched in Fig. 2 rather than the 5 SL phase that 
had been previously assumed [j, Q . Because it remains 
locally stable up to about 34.5 T (very close to the critical 
field Bc4 measured by Terada et al. [5|), the disappear- 
ance of the lOSLfTTii phase at Bed, probably occurs at a 
first-order transition between coUinear and non-coUinear 
phases. That appears to be the case for the ttii phase, 
since Bd is lower than the 7.7 T field where the SW gap 
would vanish and the tm phase would become locally 
unstable. The 10 SL TTTii phase remains locally stable 
only down to Bc2, where the frequency of a SW mode 
vanishes. 

Our work demonstrates that the stacking of the hexag- 
onal planes and the stability of a metamagnetic phase 
play crucial roles in determining the exchange and 
anisotropy parameters of a frustrated antiferromagnet. 
By constraining the fitting parameters at zero field, we 
have been able to identify the magnetic unit cell of the 
collincar metamagnetic phase in CuFe02. Constrained 
zero-field fits may prove to be a powerful technique for 
other systems as well. 
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